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ABSTRACT 

We discuss the jet kinematics of a complete flux-density-limited sample of 135 radio-loud active galactic 
nuclei (AGN) resulting from a 13 year program to investigate the structure and evolution of parsec-scale jet 
phenomena. Our analysis is based on new 2 cm Very Long Baseline Array (VLBA) images obtained between 
2002 and 2007, but includes our previously published observations made at the same wavelength, and is sup- 
plemented by VLBA archive data. In all, we have used 2424 images spanning the years 1994-2007 to study 
and determine the motions of 526 separate jet features in 127 jets. The data quality and temporal coverage (a 
median of 15 epochs per source) of this complete AGN jet sample represents a significant advance over pre- 
vious kinematics surveys. In all but five AGNs, the jets appear one-sided, most likely the result of differential 
Doppler boosting. In general the observed motions are directed along the jet ridge line, outward from the opti- 
cally thick core feature. We directly observe changes in speed and/or direction in one third of the well-sampled 
jet components in our survey. While there is some spread in the apparent speeds of separate features within an 
individual jet, the dispersion is about three times smaller than the overall dispersion of speeds among all jets. 
This supports the idea that there is a characteristic flow that describes each jet, which we have characterized by 
the fastest observed component speed. The observed maximum speed distribution is peaked at ~ 10c, with a 
tail that extends out to ~ 50c. This requires a distribution of intrinsic Lorentz factors in the parent population 
that range up to ~ 50. We also note the presence of some rare low-pattern speeds or even stationary features 
in otherwise rapidly flowing jets, that may be the result of standing re-collimation shocks, and/or a complex 
geometry and highly favorable Doppler factor. 

Subject headings: galaxies : active — galaxies : jets — radio continuum : galaxies — quasars : general — BL 
Lacertae objects : general — surveys 
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1. INTRODUCTION 

This is the sixth paper in a series in which we report results 
of the multi-epoch Very Long Baseline Array (VLBA) MO- 
JAVE program to investigate the parsec-scale jet kinematics of 
a complete flux-density-limited sample of 13 5 active galactic 
nuclei (AGN) in the northern sky at 2 cm. iLister & Homanl 
(2005) (Paper I) have described the initial source sample se- 
lection and first VLBA polarization epoch results, and in sub- 
sequent papers we have discussed the circular polarization 
characteristics (Homan & Lister! 1200 6) (Paper II), kpc-scale 
jet structure dCooper eTa l. 2007) (Paper III), and parent lumi- 
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nosity functi on (ICara & Listerl 2008)(Paper IV) of the sample. 
In Paper V (ILister et al.ll2009l) we presented a large set of 2 
cm VLBA images obtained over a 13-year period beginning 
in 1994. A primary goal of our study is a characterization of 
the long-term apparent jet motions seen in these radio-loud 
AGNs. In this paper, we discuss the kinematic properties of 
the 135 AGN jets in the flux-density-limited MOJAVE sample 
(Paper V) based on this unprecedented large VLBA data set 12 . 
We reported preliminary findings regarding a correlation be- 
tween jet speeds of the MOJAVE AGN a nd their Fermi 7- 
ray emission properties in previous papers (ILister et alj|2009t 
iKovalev et al. 2009). In subsequent papers we discuss the po- 
larization evolution, accelerations, and jet bending properties 
of the sample. 

The properties of a less well-defined sample of radio-loud 
AGNs covering the period 1994 August to 2001 March and 
through 2006 Septembe r we re discussed in earl ier papers by 
iKellermann et ali d2004l) and lCohen et alj (12007b . In these pa- 
pers we reported apparent jet velocities typically between zero 
and 15c, with a distribution extending to about 35c, and found 
that different features within individual jets appear to move 
with comparable velocity. In general, we found that lower- 
luminosity AGNs have slower jets than the strong blazars. We 
also noted the presence of apparently stationary jet features, 
as well as the presence of bends and twists where the observed 
motions do not back-trace to the base of the jet, but are instead 
aligned with the local jet direction (i.e., they are non-radial). 



12 The M OIAVE data archive 
http://www.physics.purdue.edu/MOIAVE 
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We found little or no evidence for backward motion toward 
the jet origin, such as might be expected from random pattern 
speeds. However, detailed analysis was complicated by indi- 
vidual featu res that may break up or combine with others in 
th e jet, (e.g.JKadler et al.l l2008l). 

iCohen et al.l (12007) have interpreted the observed speed dis- 
tributions in terms of the intrinsic bulk flow Lorentz factors 
and (unbeamed) synchrotron luminosities, and reported re- 
sults roughly consistent with simple relativistic beaming mod- 
els where the pattern velocity is comparable to the flow ve- 
locity. In a few cases we found clear changes in the direction 
and/or spee d of individual jet fe atures as they propagate along 
the iet (e.g jHoman et alj|2003l) . 

We chose a wavelength of 2 cm for our VLBA observa- 
tions to optimize both resolution and sensitivity, and to mini- 
mize the impact of adverse weather conditions. Complemen- 
tary radio-jet kinematic studies have been reported at shorter 
wavelen gths with better angular resolution but lower sen- 
sitivity dJorstad et al.l 1200 ll 120051) or at longer wavelengths 
with better sensitivit y to lower surface brightness features but 
poorer resolution dBritzen et alj2007ll2008HPiner et alj2007l: 
iKadler et al.ll2007l). 

Uorstad et al] (120011) used the VLBA at 7 mm and 1.3 cm 
to obtain multi-epoch observ ations of 33 7-ray blaza rs be- 
tween 1993 and 1997. Like Ke llermann et al.l d2004l) . they 
found that the jets of EGRET-detected 7-ray AGNs appear 
to be systematically faster than non 7-ray AGNs, and argued 
that superluminal jet component ejection times occur at the 
time o f 7-ray outbursts. Also using the VLBA. Uorstad et al] 
(2005) reported total and polarized 7 mm VLBA images of 
15 AGNs at 19 epochs between 1998 and 2001 and derived 
Doppler and Lorentz factors using estimates of intrinsic size 
estimated f rom the time scale o f radio outbursts. 

At 6 cm, Brit zen et alJ (120071) used an ad-hoc array and the 
VLBA to observe 293 AGNs from t he Caltech-Jodrell Flat- 
Spectrum Survey (iTaylor et al.ll 1996h . typically for 3 epochs 
spaced over the period 1990 to 2000. They reported on 
the distribution of velocities for 266 AGNs that had good 
quality multi-epoch data and also found that the most lu- 
minous sources ten d to have the fastest apparent motions 
(IBritzen et alJ l2008). 

iPiner et al.l (120071) used data from the Radio Reference 
Frame Image Database 13 obtained with the VLBA and up to 
seven other antennas at 3.4 cm to investigate jet motions in 87 
radio-bright AGNs during the period 1994-1998. Their over- 
all apparent speed distribution was peaked at low values, and 
ranged up to ~ 30c. 

The combined 2 cm VLBA Survey/MOJAVE data set pre- 
sented here and in Paper V represents a significant improve- 
ment over prior studies in terms of statistical completeness, 
temporal coverage (typically 15 epochs per source over 10 
years) and data quality (typical image resolution < 1 mas and 

rms image noise <; 0.4 Jy beam" 1 ). In this paper, we exploit 
the longer time range afforded by the newer observations as 
well as the use of extensive archival VLBA data to study the 
distribution of jet speeds in AGNs and to explore in more de- 
tail non-radial and accelerating jet flows as well as the emerg- 
ing evidence for low pattern speed features found in other- 
wise rapidly flowing jets. The overall layout is as follows. In 
Section |2] we discuss the observational set and our method 
of defining the individual jet features. In Section [3] we de- 

13 http://rorf.usno.navy.mil/RRIFD 



scribe the jet kinematic properties of the sample, and we dis- 
cuss general statistical trends in Section [4] We summarize 
our overall findings in Section [5] We adopt a cosmology with 
fl,„ = 0.27, fl A = 0.73 and H a = 71 km s" 1 Mpc" 1 . We use the 
convention S v oc v a for the spectral index. We refer to the 
radio sources throughout either by their B1950 nomenclature, 
or commonly-used aliases. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Observational Data 

The data analyzed in this study were presented in Paper V, 
and consist of a set of 2424 2-cm VLBA observations of the 
complete flux-density-limited MOJAVE sample of 135 AGNs 
spanning 1994 August 31 to 2007 September 6. Although 
the typical restoring beam of the naturally-weighted images 
has dimensions of 1 . 1 by 0.6 mas, the positions of individual 
jet features can be de term ined to much higher accuracy, as 
we describe in Section I2T2I The temporal coverage varies be- 
tween 5 and 89 VLBA epochs per source, with a median of 15 
epochs. In scheduling our observations, attempts were made 
to sample those AGNs with the most rapid jet motions more 
frequently. Some sources had considerably more archival 2 
cm VLBA data (not from our program) owing to widespread 
community interest or their popularity as calibrators. 

The MOJAVE sample is selected on the basis of compact 
(2 cm VLBA) flux density, which favors AGNs that have 
strongly beamed, fast jets viewed at small angles to the line of 
sight (i.e., blazars). The flux density limit (1.5 Jy for declina- 
tions > 0°, and 2 Jy for -20° < S < 0°) was chosen to provide 
a sufficiently large number of quasars (N = 101), BL Lac ob- 
jects (N = 22), and radio galaxies (N= 8) to investigate statis- 
tical trends in the overall sample and within the optical sub- 
classes. The redshift information on the sample is currently 
93% complete; the missing sources are mainly BL Lacs with 
featureless spectra. At the present time, only four sources 
(0446+113, 0648-165, 1213-172, and 2021+317) lack pub- 
lished optical counterparts. The full list of MOJAVE sources 
and their general properties can be found in Table 1 of Pa- 
per V. 

2.2. Gaussian Model Fitting 

The majority of the jets in our sample have relatively sim- 
ple radio morphologies that consist of a bright (presumably 
stationary) "core" and a faint jet that typically extends for a 
few milliarcseconds on one side of the core (Paper V). Given 
the modest dynamic range of a few thousand to one in our 
images, the jets can be typically modelled as a series of dis- 
tinct Gaussian features. For the purposes of tracking these 
jet features across multiple epochs, we carried out our anal- 
ysis in the interferometric (u,v) plane by fitting the bright- 
est features in each source with two-dimensional Gaussian 
components using the "modelfit" task in Difmap (Shepherd 
1997). This method offers advantages over fitting done in the 
image plane, as it uses the full resolution capability of the 
VLBA dKovalev et al.ll2.005b . However, components defined 
in this way do not necessarily correspond to real features in 
the image, but instead may merely be mathematical constructs 
needed to reprodu ce a complex brig htness distribution. In our 
previous analysis (Keller mann et alJ2004l) . we examined each 
image for continuity in position, flux density, and dimensions 
and fit Gaussian components in the image plane to determine 
component positions, thus sacrificing ef fective resolution for 
robust solutions. As discussed in Section l2.2.2l in order to in- 
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TABLE 1 
Fitted Jet Components 
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NOTE. — Columns are as follows: (1) IAU Name (B1950.0); (2) component name (zero indicates core 
component); (3) observation epoch; (4) flux density in Jy; (5) position offset from the core component 
(or map center for the core component entries); (6) position angle with respect to the core component (or 
map center for the core component entries); (7) FWHM major axis of fitted Gaussian (milliacseconds); 
(8) axial ratio of fitted Gaussian; (9) major axis position angle of fitted Gaussian. 
a Individual component epoch not used in kinematic fits. 



sure against non-physical solutions in our kinematic analysis, 
as might be introduced by (w,v) plane fitting, we have used 
only "robust" components with reliable identification across 
at least five epochs. Our fits yielded positions, sizes and flux 
densities for distinct features in the jets, which are listed in 
Table Q] 

An important constraint on our kinematic analysis is to be 
able to reliably cross-identify individual jet features from one 
epoch to the next, in order to trace their motions. A major 
difficulty arises from the fact that a model fit to an individual 
epoch is not necessarily unique. Choices can be made for the 
number of components used, and whether these are delta func- 
tions, circular Gaussians, or elliptical Gaussians. We carried 
out the majority of our fits by initially assuming an elliptical 
Gaussian for the core component, and circular components 
for the jet components. This was done to better fit unresolved 
structure very near the cores, as is sometimes the case with 
newly emerging components. We found that using circular 
Gaussians for the jet components reduced the number of free 
parameters and gave more consistent centroid positions across 
the epochs. 

After each component was fit, we inspected the residual im- 
ages for regions along the jet above ~ 5 mJy and iteratively 
added new components as necessary. If any component (in- 
cluding the core) approached zero axial ratio or an extremely 
small size during the iterations, we replaced it with a circular 
Gaussian initially comparable in size to the restoring beam, 
and then if necessary, a delta-function component. In order 
to maximize the level of continuity across epochs, we used 
the model from the previous epoch (with components restored 
back to initially circular or elliptical Gaussians) as a start- 
ing point for the next epoch. Often it was necessary to add 
or delete components as new jet features emerged and others 
faded with time below our sensitivity le vel. W e discuss the 
errors in our fitted components in Section l2.2.3l 

2.2.1. Core Component Identification 



Radio-loud blazars typically have a pc-scale morphology 
that consists of a faint jet and a bright, compact core com- 
ponent that is optically thick at centimeter wavelengths. Be- 
cause the phase self-calibration process removes any absolute 
positional information from the data, the identification of this 
component is crucial for kinematic studies, as it provides the 
(assumed) stationary reference point for all the epochs. In 
the case of the MOJAVE sample, most jets were highly core- 
dominated, which greatly simplified the identification of the 
core component. It was typically the brightest feature at the 
extreme end of a one-sided jet at most or all epoch s. Ho wever, 
in some of the two-sided jets (discussed in Section [3~6*l > and in 
a few individual sources described below, the core component 
identification was not as straightforward. 

0202+149: This source contains several jet components 
within ~ 0.5 mas of the brig ht eastern feature. H igh resolu- 
tion 7 mm VLBA images by Marscher et al. (2002) indicate a 
strong component 0.5 mas from the core, which corresponds 
to the brightest feature in our 2 cm images. Since there is 
little evidence of a counterjet in VLBA images published at 
other wavelengths, and the source is a core-dominated, mod- 
erately high-redshift quasar, we assigned the core position to 
the easternmost feature in the jet. 

0212+735: This AGN ejected a very bright moving feature 
in early 1991 that has remained significantly brighter than the 
westernmos t feature. Previous 7 mm VLBA observations by 
iListerl (120011) confirm the identification of the latter feature as 
the core. 

0316+413 (3C 84): The central feature of this two-sided 
radio galaxy contains complex, curved structure on scales 
smaller than 2 mas. There is still uncertainty as to the exact 
location of the core within th is structure, even at shorter wave- 
lengths dDhawan et al.Hl 99*81) . We have arbitrarily chosen the 
northernmost feature within the central complex as the core 
component. 

0333+321: There is significant small-scale structure near 
the westernmos t feature in the jet, b oth at 2 cm and at longer 
wavelengths dLister & Smith! 120001) . Since there is no evi- 
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dence of a counterjet in this source at 2 cm or other wave- 
lengths, we have assigned the core identification to the west- 
ernmost feature, despite the fact that is not the brightest fea- 
ture at all epochs. 

0738+313: Prior to 2001, the jet was dominated by a bright 
southern feature, which subsequently faded. We consider the 
northernmost feature to be the core of this source. 

0742+103: This is the only bona fide peaked-spectrum 
source in the flux- density-limited MOJAVE sample 
(iTorniainen et alJ 120051) . and given the known character- 
istics of this class, it is a potential candidate for two-sided 
jet structure. It is among the highest-redshift objects in the 
MOJAVE sample (z = 2.624). We have fit a slightly weaker 
component to the southeast of the brightest feature at all 
epochs, but have arbitrarily assigned the core position to the 
brightest feature in the images. We were not able to find any 
published shorter wavelength VLBI images that could be 
used to clarify the structure near the core. 

0923+392: This well-studied core-dominated blazar has a 
bright complex of eastern jet components that have been in- 
terpreted as a region where the jet b ends directly into the line 
of sight. VLBI spectral analysis bv lAlberdi et all (11997b has 
identified the (weaker) westernmost jet feature as the core. 

1226+023 (3C 273): This well-known optically low- 
polarization quasar frequently ejects very bright components 
that dominate the core emission. We have assigned the core 
position to the northeastern-most jet feature at all epochs. 

1228+126 (M87): This nearby radio galaxy has a relatively 
faint diffuse ex tension t o the east of the brightest featur e in 
the jet, which iLv et al.1 d2007|) and IKovalev et ail (12007b ar- 
gued could be a counter-jet. We have assigned the core posi- 
tion to the brightest feature, since it has the highest brightness 
temperatur e and is also the site o f weak circularly polarized 
emission (Hom an & Lister! [20061) . The kinematic fit t hat w e 
present here differs slightly from that of lKovalev et al.l (12007b . 
since we have not used the 2005 April 2 1 epoch which lacked 
the St. Croix antenna. 

1548+056: The north-south jet direction and poorer (u,v) 
coverage of this equatorial source make its core identification 
difficult. We have assigned the core to the southernmost jet 
feature. 

1958-179: There are several components within the inner- 
most 1 mas of the jet. Given the difficulty in resolving the 
emission on this small scale, we have assigned the core posi- 
tion to the brightest feature at all epochs. 

2021+614: This radio galaxy has long been considered part 
of the gigahertz -peaked spectrum class; however, owing to its 
long term variability at high frequencies and shallow spec- 
tral curvature, we have classified it as a flat-spectrum source. 
Archi val multi-frequency VLBA and space-VLBI observa- 
tions dLister et alJ^OOla ^ suggest the core is located in the 
southernmost feature, which has a flat spectral index and high 
brightness temperature. 

2128-123: This core-dominated, optically low-polarization 
quasar is too strongly v ariable for the peaked spectrum class 
(ITornia inen et al. 20 05b, despite havin g a radio spectrum that 
is slightly peaked (IKovalev et al.l 1999b . We have assigned the 
core position to the northeastern-most feature in the jet. 

2134+004: This source has a significantly curved ridgeline 
and highly complex polariza tion structure (List er & Homanl 
120051: iHoman & Listeril2006l) . Short er wavelength VLBA ob- 
servations bv lLister & Smithl (j2000) indicate that the jet core 
lies at the southeastern edge of the radio structure. 

2251+158 (3C 454.3): This powerful quasar undergoes fre- 



quent flares and component ejections that dominate the core 
emission. We have assigned the core position to the eastern- 
most jet feature at all epochs. 

2.2.2. Jet Feature Identification Across Epochs 

In most cases, the identification of jet features across 
epochs proved to be relatively straightforward, owing to the 
slow angular evolution of the source with time. We also em- 
ployed additional checks on the continuity of flux density 
and brightness temperature over time to confirm the cross- 
identifications. In some cases, however, the situation was 
much more difficult because of one or more of the following 
factors: a) insufficiently frequent temporal coverage, b) jet 
features with extremely fast angular speeds and/or flux den- 
sity decay timescales, c) a high rate of new jet feature ejection, 
d) jets with significant small scale (possibly curved) structure 
at or below the VLBA resolution level, an d e) jets con t aining 
both moving and near-stationary features. iPiner et aTl (12007b 
discuss how temporal under-sampling can often lead to com- 
ponent misidentification and incorrect estimates of apparent 
jet speeds in AGNs. 

In Table|2] we list all of the robust components that we con- 
sider to have reliable identifications across a minimum of five 
epochs. We point out that despite the fact that some features in 
the list of all fitted components (Table [TJ have been assigned 
the same ID number across epochs, only those listed in Ta- 
ble |2] have robust cross-identifications. Many of these were 
the result of comparing fits made completely independently 
by two or more authors, who were free to make their own 
choices regarding their fits and cross-identifications within the 
general guidelines described above. 

There were eight sources in the sample for which we 
were not able to measure any robust jet speeds: 0109+224, 
0235+164, 0727-115, 0742+103, 1124-186, 1324+224, 
1739+522, and 1741-038. Three of these jets were highly 
compact, with little resolved structure (see Section [3~7| i. The 
sources 0235+164, 1324+224, and 1741-038 were all ex- 
tremely compact, with no measured components at separa- 
tions larger than 0.4 mas. The components we did observe had 
considerable scatter in their fitted positions, and none were 
robust. Given that the structure in these jets was on angu- 
lar scales just at or less than the size of the restoring beam, 
we cannot quote any reliable upper limits on their apparent 
speeds. Their kinematics would be better studied with shorter 
wavelength ground and/or space VLBI, which could provide 
better angular resolution. The other five jets for which we 
were unable to measure robust speeds were: 

0109+224: This is a BL Lac object with a n unknown red- 
shift. A lower limit of z = 0.4 was claimed bv lFalomol ( 1 1996b 
on the basis of its optical host galaxy appearance. The fitted 
jet components display an unusual amount of positional jit- 
ter from epoch to epoch, which may be a result of temporally 
under-sampled fast motions. 

0727-115: There was too much positional sca tter in the 
compo nents to consider any of them as robust. IPiner et all 
(2007) analyzed 19 VLBI epochs at 4 cm and found three 
compo nents with different speed s, all with very large error 
values. Keller mann et alJ d2004l) reported a very high speed 
of 31 .2 ± 0.6 c but this was based on only three epochs at 2 
cm. 

0742+103: This source likely underwent a structural 
change in the core region in the time period 2002-2004, dur- 
ing a temporal gap in our VLBA observations. We are not 
able to measure any robust component speeds in the jet be- 
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cause of uncertainties in the core position over time. The 
parsec-scale jet structure is approximately stable since 2005, 
wh ich is consistent with the low speed (2.7 ± 0.9c) reported 
by iKellermann et al. | (120041) a t 2 cm . The three component 
speeds reported by lPiner et ail (120071) at 4 cm have large error 
values. 

1124-186: The jet structure in this low-declination source 
is very compact and extended in a southern direction that 
unfortunately coincides with the major axis of the restoring 
beam. The fitted component positions thus suffe r from high 
scatter , and we categorized none of them as robust. Pine r et all 
(120071) reported two component speeds of 1.8 ±29.5 c and 
4.0 ± 5 .4 c for this jet at 4 cm. 

1739+522: The jet is highly curved near the core, and con- 
tains considerable small scale structure that made it difficult 
to cross-identify its components across the epochs at the level 
of our angular r esolution. The single component reported by 
iPiner et all (120071) at 4 cm had an inward speed ofl6.8±7.9c. 

2.2.3. Errors on Fit Parameters 

The errors in the model fit parameters of individual Gaus- 
sian components in VLBA images are difficult to determine, 
because of non-linear dependencies on the interferometric ar- 
ray coverage, thermal noise, and the presence of other nearby 
jet features. For isolated componen ts, our previous work (see 
Appendix A of iHoman et alJ 12002) has indicated typical er- 
rors of 5% on the total intensities and positional errors of 
~ 1 /5 of the restoring beam dimension. For the brightest 
components, these errors are smaller by approximately a fac- 
tor of two. 

The vector motion (i.e., constant angu lar v elocity) and ac- 
celeration models described in Section 13.21 were fit to our 
position versus time data using a \ 2 minimization routine 
(Press et al. 1992). Given the well understood difficulties 
in obtaining robust and statistically accurate component po- 
sition uncertainties directly from the interferometric data in 
each epoch indepen dently, we followed the approach used by 
IHoman et al.l (1200 ll) of initially assuming equal weighting for 
all of our data points in a given fit, and then uniformly rescal- 
ing their weights such that \ 2 = ^points - A^ parame ters ■ In addi- 
tion to the limitation that all data points are equally weighted, 
this approach has the disadvantage that the final, minimum \ 2 
value itself is fixed to equal the number of degrees of freedom, 
and therefore it cannot be used to determine the suitability 
of the chosen motion model. However, under the assumption 
that the applied model is appropriate, this approach does yield 
good uncertainty estimates in all of the model parameters. 

We have studied the scatter of the measured component po- 
sitions about the best-fit proper motion model: this was ei- 
ther the initial vector motion fit, or the acceleration model if 
the measured acceleration was > 3a. Estimates of the typ- 
ical component position uncertainties appear in the last two 
columns of Table [2] In Figure [T] we plot the distribution of 
these x and y component position uncertainties. The indi- 
vidual distributions for x and y are so similar that we have 
simply combined the data to construct Figure [TJ We find 
the most probable component positional uncertainty is in the 
range 0.04-0.06 mas, with 68% of the components having an 
estimated positional uncertainty < 0.12 mas and 95% < 0.29 
mas. 

An altern ate means of es timating the component positional 
errors (e.g. JPiner et al.l2007l) . is to examine the differences in 
the fits to two or more individual epochs that are very closely 
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FIG. 1. — Histogram of the estimated typical position uncertainty for fit- 
ted Gaussian jet components. A total of 1052 component uncertainties are 
plotted (2 for each of the 526 robust components in our sample, one in the 
A-direction, the other in the y-direction). The large ^-position uncertainties 
for components 1 and 2 of 1228+126 (M87)(1.01 and 1.32 mas respectively) 
are not plotted. 

spaced in time, such that the jet has approximately the same 
structure in both data sets. We have done this analysis on 10 
pairs and one triplet of epochs on several MOJAVE sources 
that were within seven days of each other, and were of similar 
data quality. Based on the fastest observed proper motion in 
each source, the expected real component movement between 
the epochs in each of the pairs is less than 0.01 mas for all the 
sources except for BL Lac, which has an estimated movement 
of 0.04 mas. Therefore, the differences in the fitted compo- 
nent positions between the epochs should effectively reflect 
the positional errors of the components. 

The positional difference distribution of the individual fit- 
ted components has a narrow Gaussian central peak at zero 
and non-Gaussian tails up to roughly 0.6 mas. In terms of 
their x (R.A.) coordinates, approximately 68% of the compo- 
nents have a positional uncertainty smaller than 0.08 mas, or 
roughly 10% of the beam size. In the y (declination) direc- 
tion, approximately 95% of the components have a positional 
uncertainty smaller than 0.33 mas (roughly 30% of the beam 
size). 

We find that the positional accuracy depends on the size of 
the component and, to some degree, on its flux density (i.e., 
signal-to-noise ratio). The positional difference distribution, 
measured in units of component size convolved with the beam 
size, is very close to Gaussian. For extended components, 
a 68% (95%) confidence limit is 7% (22%) of the compo- 
nent size convolved with the beam size. For a point source 
or delta-function component with good SNR, the one-sigma 
uncertainty in position is less than one tenth the beam size. 
These error estimates are consistent with those derived from 
our multi-epoch fits as described above. 

3. PARSEC-SCALE JET KINEMATICS 

3.1. Vector Speeds 

For our initial velocity analysis of the robust features we 
have assumed non-accelerating, two-dimensional motions on 
the sky. We fit vector proper motions to the x and y positions 
of all 526 robust components in our sample, initially assuming 
simple, non-accelerating motion. For those features with at 
least ten epochs of data, we also refit the data including accel- 
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er ation terms, as desc ribed in Section [3721 and in more detail 
bv lHoman et alj (120091) . The resulting fits are velocity vectors 
relative to the component's position at a reference (middle) 
epoch listed in column 14 of Table[2] This table lists the mag- 
nitude, /t, of this vector (column 7) in microarcseconds per 
year, and its direction on the sky, (0; column 9), as well as the 
least-square fit errors on these quantities. Note that for com- 
ponents with > 3cr accelerations, the values for /i and cf> in 
Table |2] are from the fit with the acceleration terms included; 
however, as described is §3.2 these correspond closely to the 
same values obtained from the unaccelerated fit. 

In order to characterize whether the motion is directed 
purely radially toward or away from the core, we compare 
</> with the mean position angle of the component, ((i?); col- 
umn 6). We consider features in which the velocity is sig- 
nificantly offset from the purely radial (i.e., — <j>\ = 0°, or 
-0| = 180°) direction at the 3 sigma level to be "non- 
radial". These are flagged in column 10. We flag features 
as "inward" in column 9 if the velocity vector offset is sig- 
nificantly in excess of 90°. We discuss thes e rare inwardly 
moving features in more detail in Section [3~4l 

In column 13, we list the time of ejection (defined as when 
the core separation equals zero) derived from the vector fit. 
No times of ejection are listed for features in which any of the 
following apply: i) th e component shows significant acceler- 
ation (see Section [3~2l . ii) the magnitude of the vector motion 
is smaller than 3 times its associated error, or iii) the motion 
is significantly non-radial or inward. 

3.2. Accelerations 

For the features in which we have obtained at least 10 
epochs, there is sufficient positional information to attempt 
a simple tw o-dimensional accel eration fit. The method is de- 
scribed by iHoman et alj (1200 ll) . and uses the following pa- 
rameterization: 



x(0 = M*(*-*xo) + y(*-W) 2 (!) 

y(t) = (jy(t - tyo) + y (f " tmidf (2) 

Where ? m id = (?„„■„ + t max ) /2 is the middle epoch of our ob- 
servations. This epoch is chosen as the relative epoch for 
fitting accelerations, and this choice has the benefit that the 
other parameters of the fit will correspond closely to those 
from the vector motion fit (which assumes no acceleration). 
Thus, the mean angular velocities in x and y are given di- 
rectly as i(f m id) = \i x and y(f m id) = /i y respectively. The epochs 
of origin, assuming these average speeds have applied since 
the component emerged are also given directly by this fit: 

txO = *md-*ftmd)A'ftnid) and tyo = fmid-y(fmid)/.y(frnid)- 

The angular accelerations on the sky plane are given by 
X= fj, x and y = fby. The accelerations are then resolved into 
two components relative to the mean angular velocity direc- 
tion 4>: one parallel component, /iy, and one perpendicular 
component, /ij_, along 0+90°. In this way /iy directly quanti- 
fies changes in apparent speed of the component during our 
observations, and /ij_ quantifies changes in apparent direc- 
tion on the plane of the sky. In a few isolated cases where 
the angular speed of the feature is extremely slow (less than 
~ 50 /tas y" 1 , e.g., component 4 of 1253-055 (3C 279), there 
may not be a well-defined velocity vector direction. In such 
cases, the individual perpendicular and parallel components 



of the acceleration are not necessarily meaningful. 

We note that our approach to fitting for angular accelera- 
tions assumes that the same acceleration has applied through- 
out our observing period. Other models, such as impulsive 
accelerations inducing sudden changes in speed and/or direc- 
tion are certainly possible and may provide a better fit in some 
cases, such as component 1 in 1253-055 (3C279) which was 
re ported to make such a sudden change in speed and direction 
bv lHoman et all (I2003T) . 

We list the parallel and perpendicular components of the 
fitted accelerations in columns 11 and 12, respectively, of Ta- 
ble|2] We have flagged in column 7 any features which have a 
parallel or perpendicular acceleration exceeding the 3 sigma 
level. Of the 311 features we examined, 109 (35%) fall into 
this category. The angular speed (p) listed in Table|2]for these 
features reflects the value obtained from the acceleration fit, 
as opposed to the vector one. It is important to note that there 
are likely other jet features in Table|2]that are in fact changing 
speed and/or direction, but we cannot discern them as such 
because of insufficient temporal coverage. 

In Figure |2] we plot our fits to each source in a core separa- 
tion distance versus time diagram. Each robust jet component 
is indicated by the colored points, which are joined by either 
a solid line (non-accelerating fit) or a dotted line (accelerat- 
ing fit). Thin line-widths are used in cases where the compo- 
nent motion is significantly non-radial. Unfilled colored cir- 
cles indicate individual data points that were not used in the 
kinematic fits, and unfilled black circles indicate non-robust 
components. 

We plot the motion fits and sky positions of the individ- 
ual robust components for selected AGN in Figures [3] and [4] 
Extended versions of Figures [2] [3] and [4] containing plots of 
all of the robust jet components in the MOJAVE sample, are 
available online. The left hand panels of Figures[3]and|4]show 
a 2 cm VLBA contour image of the source from Paper V at 
the middle epoch listed in Table |2] The orange box delimits 
the zoomed region displayed in the right hand panels. The 
component's position at the middle epoch is indicated by the 
orange cross-hairs. The dotted line connects the component's 
middle epoch position with that of the core feature. The po- 
sitions at other epochs are plotted with unfilled blue circles. 
The red solid line indicates the vector (or accelerating) fit (see 
Table [2) to the component positions. The red arrows in the 
right hand panels indicate the direction of motion, and the 
gray dashed circles/ellipses represent the FWHM sizes of the 
individual fitted Gaussian components. 

3.3. Non-radial Motions 

In our previous study of AGN jet kinematics in the 2 cm 
Survey (Keller mann et al.l 120041) . we identified many cases 
where jet features were moving in directions that did not ex- 
trapolate back to their presumed origin at the base of the 
jet. These non-radial components provided convincing evi- 
dence that jet features change direction, and ruled out sim- 
ple ballistic "cannon ball" models. In our study of the MO- 
JAVE sample, we find that 166 out of the 526 robust compo- 
nen ts (3 2 %) satisfy our definition of non-radial motion (Sec- 
tion l3.U . A comparable fra ctio n was found i n sma ller samples 
by Kellerm ann et al.l (120041) and lPiner et al.l (|2007), indicating 
that non-radial motion is a common feature of the jet flow in 
blazars. We note that three of these non-radial components 
might be considered "inward", i.e., their direction of motion 
is greater than 90° from the local jet direction at the 3cr level. 
Five additional components that would qualify as inward do 
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FIG. 2. — Plot of angular separation from core versus epoch for Gaussian jet components. The B1950 source name is given at the top left of each panel. 
Color symbols indicate robust components for which kinematic fits were obtained (dotted and solid lines). The solid lines indicate vector motion fits to the data 
points assuming no acceleration, while the dotted lines indicate accelerated motion fits. Thick lines are used for components whose fitted motion is along a 
radial direction from the core, while the thin lines indicate non-radial motions. Unfilled colored circles indicate individual data points that were not used in the 
kinematic fits, and unfilled black circles indicate non-robust components. The component identification number is indicated next to the last epoch of each robust 
component. (This is a figure stub; an extended version is available online.) 



not qualify as non-radial since their | (i)) - (j>\ values are less 
than 3a from 18 de grees. We discuss these features in more 
detail in Section [3~4l 

The high incidence of non-radial motion among the com- 
ponents in our sample indicates that changes in jet motion are 
common. A mo re detailed analysis o f a subset of these data, to 
be presented by Ho man et alj {2009), indicates that in the vast 
majority of cases, the non-radial motions are indeed caused 
by accelerations in a direction perpendicular to the velocity 
vector. 

3.4. Apparent Inward Motions 

For the great majority of our reported jet features, we find 
that their motion is either directed outward from the core, or 
is nearly stationary. Fewer than ten of the 526 robust features 
we have studied show evidence for inward motion, i.e., where 
the projected separation between the jet feature and the core 
decreases with time, and most of these are only marginally 
significant. The illusion of apparent inward motion can be 
explained in terms of one or more of the following scenarios: 

(i) If a newly emerging jet feature and the core are not re- 
solved by our beam, this would cause a temporary shift in 



the measured position of the core centroid and a correspond- 
ing apparent decrease with time in the apparent separation of 
the co re and other more distant jet features (e.g., Gom ez et"aT] 
fl995l) . 

(ii) It is also possible that the true core is not seen, possi- 
bly due to absorption, and that the only observed features are 
parts of a jet w hich are moving with different velocities (e.g., 
Lobanov 1998). If the furthest component is moving with a 
slower velocity than the one closest to the obscured core, then 
the two components will appear to be approaching each other. 

(iii) An apparent decrease in component separation from 
the core could also arise due to component motion away from 
the core but along a highly curved jet. If the jet bends back and 
across the line of sight, the projected separation from the core 
can appear to be decreasing, even though the actual distance 
from the core is increasing with time (e.g., Marsch er et alj 
IT991 . 

(iv) Some observed jet features may merely reflect patterns 
in the flow, of which so me might even be moving backward. 
Ilstomin & Parievl JT996) have discussed an electron-positron 
jet model where an observer located close to the d irectio n of 
motion will see backward-moving knots. iGongl d2008l) has 
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FIG. 3. — Vector motion fits and sky position plots of individual robust jet features in MOJAVE AGN. Positions are relative to the core position. The left hand 
panels show a 2 cm VLBA contour image of the source from Paper V at the middle epoch listed in TableffJ The orange box delimits the zoomed region displayed 
in the right hand panels. The component's position at the middle epoch is indicated by the orange cross-hairs. The dotted line connects the component with the 
core feature and is plotted with the mean position angle {1?) (Tableff). The position at the middle epoch is shown by a filled violet circle while other epochs 
are plotted with unfilled blue circles. The red solid line indicates the vector (or accelerating) fit (see Tableff) to the component positions. The red arrows in the 
right hand panels indicate the direction of motion, and the gray dashed circles/ellipses represent the FWHM sizes of the individual fitted Gaussian components. 
Displayed from top to bottom in the figure are component ID = 14 in 0415+379 (3C 111), ID = 3 in 0738+313, and ID=6 in 0738+313. 
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FIG. 4.— Continued. Displayed from top to bottom in the figure are component ID = 4 in 0829+046, ID=5 in 1253-055 (3C 279), and ID=2 in 2351+456. (This 
is a figure stub; an extended version is available online.) 
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proposed a model of a precessing jet which illuminates an am- 
bient screen, so that an observer oriented near the plane of the 
screen can see projected motion in any direction even though 
the jet itself is always flowing outward. We note, however, 
that any model involving patterns must account for the the 
very large numbers of observed outward motions, as opposed 
to inward ones. 

(v) Internal changes in the brightness distribution of a large, 
diffuse feature may cause the centroid of the fitted Gaussian 
to wander, possibly mimicking inward motion. 

We discuss several individual cases for inward motion in 
our sample below: 

1458+718: This compact steep-spectrum quasar has a very 
long jet with a bright complex region located about 25 mas 
from the nucleus as well as a fainter region about 40 mas 
downstream. Close to the nucleus, as well as for the more 
distant feature, the motion is in the outward direction with 
Amp ~ 7. It is more difficult to accurately determine the mo- 
tion of the more diffuse and complex feature located about 
25 mas downstream, but there may be a slight gross inward 
motion with /3 app ~ 2. 

1803+784: This is a classical core-jet source. The most 
distant, somewhat extended, feature is moving away from the 
core with a velocity of 9.0 ± 2.5c, and two intermediate fea- 
tures appear almost stationary over the past ten years or more. 
The two innermost features each appeared to be moving away 
from the core with a velocity of ~ 1 .5c until the year 2000, 
when they both appeared to slow down. Since 2004, the dis- 
tance between of each of them and the core has decreased. 
The very similar apparent motion of these two inner features 
is characteristic of an apparent change in separation due to an 
emerging, yet unresolved new component. However, no such 
feature is evident in VLBA data obtained as of March 2009. 

2005+403: There is a pronounced bend in the jet located 
about 2 mas from the nucleus. Close to the nucleus, there is a 
bright feature which has moved 0. 1 mas closer to the nucleus 
over the three year period that it has been observed, but this 
result depends primarily on the first and last observations in 
2004 and 2007 respectively. Between mid 2005 and mid 2006 
there was little or no motion. If real, the apparent inward 
motion of this feature is most likely due to a new emerging 
component which is unresolved from the core. Beyond about 
1 mas, the flow is in the general direction of the jet ridge line. 

2021+614: This jet has a slow-moving feature (ID = 2) lo- 
cated approximately 3 mas downstream from the core that is 
classified as inward according to our formal definition. The 
fitted angular speed is extremely slow (1 1 ± 3 /xas y" 1 ), and 
although there is a general downward trend in its separation 
versus time plot, the component oscillates about its mean po- 
sition over time. 

2200+420: The innermost component (ID=7) of this source 
changes position angle from -160°to ~ -170° from 1999 to 
2006, while maintaining a relatively stable separation of 0.3 
mas. Since this is at the extreme limit of our resolution, and 
the component has a low speed (Section [3.5t , the vector mo- 
tion direction is not necessarily well-defined. We thus cannot 
conclude for certain if the motion is truly inward. 

2201+171: Close to the core, the flow is along the direction 
of the jet, but 1.4 arcseconds downstream, a bright feature 
appears to flow in the reverse direction, possibly because of 
confusion with a newly emerging component. 

2230+1 14 (CTA 102): This AGN has an elongated jet with 
several pronounced twists and bends. Along most of the jet, 



the flow appears to be close to the jet direction, with apparent 
velocities between c and 16c. A single component at ~ 6 mas 
from the core is moving radially inward at 0.97 ± 0.22c. The 
contour maps indicate an apparently sharp bend of approxi- 
mately 50° in the jet at this location. 

2351+456: As shown in Figure setQ] 2351+456 has a well 
defined jet close to the core which has a sharp nearly right- 
angle bend about 7 mas downstream, beyond which it opens 
into a very diffuse structure characteristic of a direct head-on 
jet. Close to the core and near the bend the flow follows the 
ridge line with an apparent /3 app of 10 to 15. However a bright 
feature located about 4 mas downstream is moving toward the 
core with /3 app = 28.3 ± 1.3. Over the 6 year period that we 
have observed this feature it has moved a total of about 2 mas, 
so the apparent motion cannot be explained by the emergence 
of a new feature unresolved from the core. Rather, we suggest 
that the jet of 235 1+456 is viewed nearly end-on, and that the 
region near 4 mas is where the jet bends back across the line 
of sight. VLBA data from the MOJAVE program obtained 
after 2007 indicate that the component at 4 mas has stopped 
moving inward, supporting this interpretation. 

3.5. Dispersion of Apparent Speeds Within Individual Jets 

Past multi-epoch VLBI studies of individual radio-loud 
AGNs have identified instances where the jet displ ays features 
that move with a range of apparent speeds (e.g., iPiner et all 
120071 lJorstad^raT]l200l IKellermann et al.1 120041) . Such in- 
stances may suggest that the motions of some jet features 
are not related to the flow, or that the jet is experiencing a 
change in geometry over time due to precession. Some jet 
features have also been seen to either remain virtually sta- 
tionary over time, or move at significantl y slower apparent 
speed s than other features in the same jet (IKellermann et al.l 
12004 iBritzen et all 120051: Uorstad et all 120051) . These low- 
pattern speed (LPS) features have sometimes been explained 
as discrete locations in a bent jet that are at an extremely 
low viewing angle b ecause of the specific geometry involved 
(lAlberdi et al.ll2.000T) . In the specific case of a flow moving 
relativistically along a helical path, the changing viewing an- 
gle is predicted to give rise to Doppler boosted regions that 
appear nearly stationary in VLBI images. Alternatively, it is 
possible that some LPS features represent standing shocks in 
the flow. These can take the form of re-collimation shocks in 
an initially over-pressurized outflow, and have been success- 
fully reproduced in numerical simulations o f AGN jets (e.g., 
iGomez et all 119951: iPerucho & Martil 12007b. Some studies 
(e.g.. iMandal & Chakrabartil 120081: iBecker etal.1 120081) have 
suggested that standing shocks within the unresolved cores of 
AGN jets may play a major role in accelerating particles near 
the base of the jet, and could be respons ible for the persis- 
tent high levels of polarizat ion in blazars (jD'Arcangelo et al. 
l2007llMarscher et al.ll2008l) . 

We have examined the distribution of apparent speeds 
within each jet in the MOJAVE sample, and find a median 
rms dispersion of RMSj e t = 0.065 mas y" 1 , or 2.6c. We also 
calculated the median speed value of each MOJAVE jet, and 
find that the distribution of this value for the entire sample has 
an rms dispersion of RMS samp ie = 0.29 mas y" 1 , or 7.3c. The 
RMS samp i e and RMSj e t values differ by mo re than a factor of 
three. This supports the previous finding of IKellermann et all 
(2004) that most jets generally eject features with a character- 
istic speed. We found even larger differences between these 
two RMS values in the cases where i) the LPS features were 
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FIG. 5. — Distribution of projected distance from the core for low pattern 
speed jet features in the MOJAVE sample. 

dropped from the analysis, or ii) only quasars were consid- 
ered, or iii) the maximum jet speeds were used to calculate 

RMS samp l e . 

We have also used the regular time coverage and complete- 
ness of the MOJAVE sample to investigate how prevalent LPS 
features are in blazars, and whether they are more prone to oc- 
cur in certain types of AGN jets. We examined all 94 robust 
features in our sample with radial speeds less than 50 /ias y" 1 
to see if their motions were non-accelerating and significantly 
slower than other features in the same jet. We were able to 
identify only 31 such instances, which are indicated with a 
flag in column 5 of Table [2] These LPS features have a bi- 
modal distribution, with the majority lying within projected 
distances of 0.8 mas (~6 parsecs) from the core, while only a 
handful are found further downstream (Figure [5]). 

Since the close-in LPS features are located approximately 
within one beam-width of the core, we have investigated the 
possibility that they may be artifacts arising from the mod- 
elfitting process. We model fitted the sources using different 
models for the core, including circular and elliptical Gaus- 
sians, as well as delta-functions, and were not able to elimi- 
nate the close-in features. They generally are present in all the 
epochs on a particular source, except on occasions when they 
are blended with new components that have emerged from 
the core. After such components move downstream, the sta- 
tionary feature then reappears, indicating that it is a persistent 
feature of the flow. We have a lso compared our m odel fits 
to five jets in common with the iJorstad et al.1 (120051) sample, 
who model fitted several closely-spaced VLBA epochs at 7 
mm with three times better angular resolution. In each case, 
they identified a stationary component at roughly the same 
core distance as in our fits (taking into account a typical core 
shift b e tween these wavel engths of ~ 0.1 mas. iKovalev et all 
20081). IJorstad et alJ (12001 found in some cases additional 
stationary features even closer to the core (r ~ 0.1 mas), 
which are below the resolution level of our 2 cm images. De- 
spite our limited sensitivity to these extreme close-in features, 
we conclude that LPS features are very rarely found in blazar 
jets at projected distances greater than 6 pc downstream of 
the core. Given the small viewing angles of the MOJAVE 
jets, this translates into an approximate de-projected distance 
of50pc. 

Given the caveats of our limited angular resolution, we find 
that the incidence rate of LPS features is much higher for the 



BL Lacs in the sample (10 of 22), as compared to the quasars 
(15 of 86). Although the latter generally have higher redshifts 
than the BL Lacs, we find that for the whole MOJAVE sam- 
ple, K-S tests indicate no significant differences in the red- 
shift or 2 cm VLBA luminosity distributions of the AGN jets 
with LPS features and those without. It is therefore possible 
that the conditions in BL Lac jets are more conducive to the 
formation of strong, stable re-colli mation shocks. Numerica l 
hydrodynamicaljet simulations by Go mez et al.l d 1 9951 1 1 9971) 
have shown that the strength and separation of jet knots as- 
sociated with re-collimation shocks is dependent on the jet 
opening angle and Mach number, as well as the external den- 
sity gradient. 

3.6. Two-sided Jets 

Although it is thought that quasars and radio galaxies have 
intrinsically two-sided jets which feed the distant radio lobes, 
the great majority of the sources in the MOJAVE sample have 
one-sided jets, which is not unexpected considering the ex- 
pected large differential Doppler boosting of the approaching 
and receding jets. A few sources, however, show two-sided 
jets; these are mostly low luminosity radio galaxies, and it is 
tempting to consider that these sources lie close to the plane 
of the sky so that there is littl e differential Doppler boosting. 
However, ICohen et alJ d2007l) have shown that the observed 
number of slowly moving low luminosity two-sided jets is 
considerably greater than would be expected if the distribu- 
tion were determined by orientation alone, so that the two- 
sided jets may reflect a population intrinsically different from 
the one-sided population. 

One challenge in determining independent speeds in a two- 
sided jet is determining the unambiguous location of the com- 
mon, presumably stationary, central core presumed to lie at 
the base of each jet. We have been able to obtain robust speed 
measurements of both jet and counterjet components for five 
two-sided jets in our sample. Stacked images over all epochs 
for these two-sided sources are shown in Paper V and are dis- 
cussed below. 

0238-084 (NGC 1052): This is a bright nearby el- 
liptical g alaxy that harbors a well-studied low luminos - 
ity AGN. IVermeulen et ail (120031) : iKellermann et alJ (120041) : 
Kadle r et alJ (12004b have previously discussed the observed 
outflow in two nearly symmetric jets containing multiple fea- 
tures which propagate away from a central region with rel- 
atively slow subluminal velocities. VLA observations show 
the jets extending symmetrically o ut to 1 .5 kpc on either side 
(ICooper et alJl2007tlWrobe]||19 84). The central region is sur- 
rounded by a rich molecular (OH and H2O), atomic, and 
ionized gas cloud. HI observations indicate the presence of 
three velocity streams, and multi-frequency VLBA observa- 
tions show patchy free-free absorption, which may come from 
the torus surrounding the central black hole (IVermeulen et alJ 
12003b iKadler et all 12004). The analysis of the jet motions 
is complicated by two factors. First, the presumed base of 
the jets is obscured by free-free absorption, especially on 
the receding side, and second, the simultaneous presence of 
many features in each jet exacerbates the problem of cross- 
identifying components over different epochs. We have esti- 
mated the location of a virtual center by assuming it is located 
half way between pairs of symmetrically moving components 
in opposing jets. Figure set[3]shows the location of each fea- 
ture that we were able to identify as a function of time. Over 
the past decade, a new feature is ejected at an average rate of 
about once per year, each moving with an apparent velocity 
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in the range of 0.1c to 0.25c. 

0316+413 (3C 84): NGC 1275 is a prominent radio 
galaxy that contains a bright complex central core, a south- 
ern cocoon-like expanding cloud which is connected to a cen- 
tral feature by a thin faint jet, as well as an inverted spect rum 
counter feature dDhawan et al.lll998t I Walker et all 120001) . It 
was detected in gamma-rays during the first 4 months of 
the Fermi observatory, during which time 2 cm MOJAVE 
VLBA obs ervations indicated a significant brightening of the 
core region lAbdo et alJ J2009). We have determined the mo- 
tion in the southern lobe to be /i = 0.27 ± 0.05 mas y" 1 (v = 
0.31 ± 0.06 c) while the northern lobe appears to move more 
slowly with fi = 0.07 ± 0.02 mas y" 1 (0.08 ± 0.03 c). Assum- 
ing the jets are intrinsically symmetric, these observed speeds 
correspond to a Lorentz factor of 0.6 with the motion oriented 
at about 1 1 degrees to the line of sight. However, the central 
component itself has a compl ex morphology, wi th observed 
motions up to 0.2c or more (Dhawan et al. 1998), which in- 
troduces some additional uncertainty in the derived motions 
of the two outer features. We have taken the origin of motion 
to be at the northern end of the central feature, as suggested 
by the 7 mm VLBA i mages and correspo nding motions within 
the central feature (Dhaw an et al.| [l998). 

1228+126 (Virgo A): M87 harbors a prominent jet which 
has been imaged at radio, optical, and X-ray wavelengths. 
We have taken the core to be coincident with the brightest 
radio component, which is sup ported by the circular polariza- 
tion reported at this position (Homan & Lister 2006). We are 
able to identify six features along the jet and one counterjet 
feature which have existed continuously over the decade-long 
interv al covered by our o bservations. As we reported previ- 
ously dKovalev et alJl2007h . we find no evidence for any mo- 
tion faster than a few percent of light speed, with the fastest 
feature measured at 0.05c ± 0.02c. The counterjet may be 
moving at about 0.01c, but this is not well determined be- 
cause of its complex structure. Faster an gular motions corre- 
sponding to ~ 2c have been reported by lAcciari et alJ (l2009t) 
in closely temporally-spaced 7 mm VLBA observations, but 
these are still not sufficient to completely account for the large 
jet-to-counterjet ratio in regions 3 mas downstream of the core 
by Doppler boosting alone. Therefore, the bulk plasma flow 
speed must be much grea ter than the observed pattern mo- 
tions dKovalev et alJ l2007). and there is the additional possi- 
bility that the M87 jet is intrinsically asymmetric. We note 
that the extreme proximity and low luminosity of M87 set it 
apart from the other AGN in our sample, which are typically 
powerful FR-II class radio sources. 

1413+135: This is a complex, highly unusual two-sided 
BL Lac object located within an apparent spiral host galaxy 
dVeron-Cettv & Veronl I2006T) . Like NGC 1052, the AGN 
is s urrounded by a cloud of molecular, atomic, and ionized 
gas ((Permian et al. 2002), but 1413+135 is a thousand times 
more luminous than the AGN in NGC 1052. We have iden- 
tified seven robust jet features in 1413+135, two of which 
are located toward the north east in an apparent counterjet. 
Both counterjet features are weak, and so it is difficult to 
determine their motions with high precision. Based on ten 
measurements over a seven year period, we find a best-fit 
velocity of the more distant counterjet feature (ID = 1) of 
1. 34 ± 0.56c, but we cannot exclude the possibility that this 
feature is stationary. The other counterjet feature (ID=2) is 
slow (0.14 ±0.28) and is also possibly stationary. By com- 
parison, in the main jet, there are several bright, robust fea- 



tures with apparent speeds ranging from near stationary for 
the strongest feature to 1.8 ±0.2c for the westernmost com- 
ponent (ID = 3). 

1957+405 (Cygnus A): Cygnus A is a classical (archetypi- 
cal) FR II radio galaxy. The central compact radio source lies 
at the origin of the well known prominent kiloparsec jet struc- 
ture. It has a luminosity comparable to other powerful AGNs, 
but like the M87 AGN, the Cygnus A AGN is many orders 
of magnitude less luminous than the extended jets and lobes. 
Also, as in M87, the jets are asymmetric, with the western jet 
being brighter on scales of few parsecs to a hundred kilopar- 
secs. Interestingly, the ejection dates of the counterjet com- 
ponents match up well with their counterparts in the main jet, 
indicating near-simultaneous ejections. On the bright, pre- 
sumably approaching side, the jet appears to be moving with 
a velocity close to 0.2c and on the faint receding side, 0.1c. 
Assuming that the two sides are intrinsically similar, the in- 
trinsic velocity is about 20 percent of the speed of light, with 
the jets oriented at an angle near 70 d egrees to th e line o f sight. 
Multi-wavelength VLBA studies bv lBach et ail d2005) report 
evidence for an absorbing feature in front of the receding jet. 

All of the apparent two-sided jets in our flux-density- 
limited sample are found in relatively low luminosity and rel- 
atively nearby radio galaxies. We have not found any sym- 
metric parsec scale jets in any quasars, although 101 out of 
the 135 sources in our sample are classified as quasars. These 
two-sided jets provide additional evidence that the apparent 
motions do indeed reflect the underlying jet flow and are not 
simply a (shock) pattern motion which propagates within a 
relatively stable or slowly moving plasma. In the case of ran- 
dom pattern speeds, we should not expect motions in the jet 
and counterjet to be correlated, whereas if the symmetric out- 
flow is due to a common explosive origin the kinematics of 
each side might show similar properties. In at least NGC 1052 
and Cygnus A we can make a one-to-one association between 
features on the two sides, supporting the idea that they re- 
flect the real jet flow, although in M87 no such association 
can be made. We note that in all of these sources, the intrinsic 
flow speeds are only mildly relativistic, and the spatial scales 
are much smaller than the more distant objects, thus we can- 
not make any inferences from this argument about the nature 
of the apparent superluminal speeds that are observed in the 
more powerful quasars. 

3.7. Highly Compact Jets 

We have also identified three extremely core-dominated 
sources (0235+164, 1324+224, 1741-038) which had very lit- 
tle discernible mas-scale jet structure at any epoch. This could 
either be because their jets are atypically weak, or they may 
possess a highly beamed core which swamps the jet emission 
in our limited dynamic range images. It is also possible that 
they may be experiencing a period of inactivity in which new 
jet components are not being ejected. The latter was the case 
for 1308+326, which had very little discernible jet structure 
prior to 1997, and now has a pronounced mas-scale jet. 

With the exception of 0235+164, none of the three sources 
show any extended radio structure in deep VLA image s 
(lUlvestad et alJ[l98ll: IPerlev et alJ[l98l ICooper et alJl2007h . 
They are therefore well-suited as calibration sources for flux 
density and absolute polarization position angle scaling when 
simultaneous single-dish/VLA and VLBA data are available. 

3.8. Stable Jets 
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log angular speed (mas y" 1 ) 

FIG. 6. — Angular speed distribution for all 526 robust jet features in the 
MOJAVE AGN sample. The shaded histogram represents the distribution of 
the fastest robust feature in each jet. 

VLBI sources that maintain a stable mas-scale structure 
for long periods with well-defined centroids are very useful 
in geodetic appl ications such as the definition of the ICRF 
( Feiss el-Vernierll2()03l) . Our long-term VLBA study has indi- 
cated, however, that such sources are exceedingly rare. We 
have identified only two sources in the MOJAVE sample 
whose components have maximum apparent speeds smaller 
than 10 pas y" 1 . Both of these are at hi gh re dshift. Unlike 
the compact sources discussed in Section 13.71 which may be 
in temporarily quiescent states, the decade-long stability of 
these sources suggests that they have low intrinsic expansion 
velocities and/or high jet viewing angles. However, we cannot 
exclude the possibility that faster changes of the jet structure 
of these quasars might occur in the future. 

0552+398 (DA 193): This AGN (z = 2.363) has been fre- 
quently classified as a gigahertz-peaked spectrum source in 
the literature. However, it displays consid erable flux den- 
sity v ariability above its turnover frequency (Torniain en et alJ 
2005), and has persistent linear polarization on mas scales that 
is uncharacteristic of the GPS class. We were able to fit all of 
the 2 cm MOJAVE VLBA epochs between 1997 and 2007 
with a two component Gaussian model, in which the sepa- 
ration is barely increasing at a rate of 4 ± 1 /ias y" 1 (0.36 ± 
0.09c). 

0642+449: This high redshift (z = 3.396) quasar has also 
been p reviously classified as a GPS, yet iTorniainen et alJ 
(2005) found its flux density variability to be too high for this 
class. During our 12 year monitoring interval we have found a 
single robust jet component whose core separation is increas- 
ing at only 7±1 ^asy -1 (0.76 ± 0.11c). 

4. OVERALL KINEMATICS OF THE MOJAVE JET 
SAMPLE 

4. 1 . Distribution of Angular Jet Speeds 

The robust jet features in our sample span a wide range of 
apparent speed, measured in both angular units and in units of 
the speed of light. The distribution of angular speeds (Fig.|6]l 
provides useful information for the pl anning of kinemati c 
monitoring programs (e.g., TANAMI: iKadler et all 120071) . 
since it indicates how often bright blazar jets should typically 
be observed to accurately track their motions. The distribu- 
tion is roughly log-normal, with a median of 0.22 mas y -1 . 
The feature with the largest measured proper motion is in 



the jet of the nearby radio galaxy 0430+052 (3C 120) and 
is moving at 2.46 ±0.074 mas y" 1 . In the MOJAVE survey, 
we have attempted to observe sources such that individual jet 
features move less than a beam width between observations. 
This facilitates the identification of features across successive 
epochs. Cadences range from one observation per month for 
the fastest sources to once every two to three years for the 
slowest sources, with a median of once every 1 1 months (Pa- 
per V). 

4.2. Maximum Jet Speeds 

A primary goal of the MOJAVE survey is to measure the 
distribution of jet speeds in our flux-density-limited AGN 
sample, in order to investigate the overall speed distribu- 
tion and beaming properties of the blazar parent population. 
Such information is needed for studies of AGN luminosit y 
functions and blazar demographics (e.g., Cara & Lister 2008). 
The issue is complicated somewhat by the possible presence 
of multiple pattern speeds in an individual jet, some of which 
may be moving at a different speed than the bulk flow (Sec- 
tion [33}. We are therefore interested in the fastest speed in 
the jet measured over our approximately decade-long moni- 
toring period. For the purposes of constructing a maximum 
speed statistic for each source, we considered all robust com- 
ponents in Table [2] with angular speed measurements above 
the 3cr level. This method addressed some isolated cases such 
as2201+171, where a diffuse component with large positional 
uncertainties gave a large (and poorly determined) speed mea- 
surement that dominated the statistic for the source. In the 
case of 0422+004, which had only a single robust component, 
we adopted its speed as the maximum speed, despite its mea- 
sured value being slightly below 3 a. 

4.2.1. Distribution of Maximum Jet Speeds 

In Figure [7] we plot the distribution of jet speeds for all ro- 
bust features in the sample. The shaded histogram represents 
the distribution of maximum speed. The fall-off towards high 
values is consistent with previous large VLBI surveys car- 
ried out at longe r observing wavelengths (Britze n et alj|2007t 
iPiner et all 12007b . However, a major difference in the MO- 
JAVE distribution is that it is peaked at a moderately high 
value (/3 a pp — 10), whereas the earlier surveys were peaked 
at values < 4c. The selection criteria of all three surveys are 
similar, and contain many overlapping objects. The MOJAVE 
survey does have higher angular resolution and excellent tem- 
poral coverage, however, which may facilitate the identifi- 
cation of fast-moving components that would otherwise be 
missed. Indeed, our fastest identified jet component (compo- 
nent ID=3 in the quasar 0805-077 14 ) at 50.6 ±2.1 cis signifi- 
cantly faster than any previously identified in cm-wavelength 
VLBA surveys, and is comparable to some of the fastest com- 
ponents identified in a small sample of 15 blazar jets obs erved 
monthly with the VLBA at 7 mm by lJorstad eTaTI (120051) . The 
median measured speed in the latter study was also ~ 10c, 
which suggests that high angular resolution is essential for 
determining accurate jet speeds in blazar proper motion sur- 
veys. Data from the Fermi observatory have shown t hat the 
extremely fast MOJAVE blazars, including 0805-077 dCiprinil 

14 The r e dshift of 0805—077 (z = 1.817) was originally measured by 
IWhite et"aT] (1988), on the basis of several broad emission lines in the op- 
tical spectrum, and was later confirmed at the Nordic Optical Telescope by T. 
Pursimo (2009, private communication). 
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FIG. 7. — Distribution of apparent speed for 502 robust jet features in MO- 
JAVE AGN with measured redshifts. The shaded histogram represents the 
distribution of the fastest robust feature in each jet. 



2009ab), show a higher tendency to emit high energy 7- 
rays, likely as a result of their high Do ppler boosting factors 

dLister et alj200^lKovdev"etd]|2009b. 

Other studi es dVermeulen" & Cohen [1991 

iLister & Marscherl 1 1997b have investigated how the dis- 
tribution of fastest apparent speeds in an AGN sample 
selected on the basis of beamed jet emission reflects the 
overall distribution of Lorentz factors (T) in the parent 
population. First, the maximum observed apparent speed sets 
an approximate value for r max , since T > (/3 app + l) 0,5 , and in 
a large parent population, some of the jets with T = r max will 
have their observed flux densities Doppler boosted by very 
large factors by virtue of their near end-on orientations. For 
reasonable choices of the parent luminosity function, most 
of these sources will be selected in a flux-density-limited 
survey, regardless of their redshift. Second, the shape of 
the distribution in Figure [7] rules out a single-valued T 
distribution, since in that case the apparent speed distribution 
would be expected to peak at ar , p g T, and drop off s harply 
for smaller values of /3 app dVermeulen & Cohenlll994l) . The 
tapering of the speed distribution at higher values strongly 
suggests a power-law Lorentz factor distribution, which 
we use in M onte Carlo beaming simulations described in 
Section l4T2~3l 

4.2.2. Maximum Speed Versus Redshift 

In Figure [8] we plot maximum speed against redshift for 
the full sample. The data points are approximately evenly 
scattered for redshifts greater than ~ 0.6, whereas at lower 
redshifts, there is an apparent upper envelope, indicating that 
the maximum superluminal speed increases with redshift. 
The dashed curve represents a constant angular velocity of 
2.5 mas y -1 , which corresponds to the largest angular speed 
measured in the survey. Since this curve lies well above the 



data points past z = 0.1, we conclude that the upper envelope 
is not the result of survey measurement bias. 

The envelope instead likely arises because of the com- 
bined effects of relativistic beaming and the sample's lower 
flux density cutoff. In the MOJAVE survey, the minimum 
detectable luminosity rises sharply wit h redshift, creating a 
classical Malmquist bias (see Fig. 1 of Cara & Lister 2008). 
The high redshift sources therefore have higher apparent lu- 
minosities, which they achieve primarily via Do ppler boost- 
ing. Our Monte Carlo simulations (Section r4.2.3t suggest that 
these jets have some of the highest Doppler boosting fac- 
tors in the parent population (i.e., ~ 2r max ), and therefore 
can have apparent speeds up to /3 app = r max . At lower red- 
shifts (below z — 0.6), beaming is not as important a fac- 
tor in determining whether a jet makes it into the sample, 
and the available co-moving volume is much smaller. This 
makes it far less probable that highly superluminal jets will 
be seen at low redshifts. The transition at which this oc- 
curs corresponds to a redshift z where the survey luminos- 
ity limit roughly coincid es with a flatt ening in the beamed 
luminosity function (LF). iListerl ((2003) showed that the lat- 
ter occurs at Li<Sax> where L\ is the low-luminosity cutoff 
of the intrinsic LF, <5 max is the maximum Doppler boosting 
factor, and p is the boosting index. For a flat spectrum jet 
with typical boosting index of p = 2, <5 max = 2/3 app max = 100, 
andLi = 1.2 x 10 23 W Hz" 1 , which corresponds to the faintest 
source in the survey (1228+126=M87), the transition occurs 
at z — 0.55, in approximate agreement with Figure[8] 

4.2.3. Maximum Speed Versus Apparent Radio Luminosity 

In an earlier p aper using preliminary MOJAVE data, 
ICohen et al.l d2007l) reported an upper envelope in a plot of 
maximum jet speed versus VLBA luminosity. In Figure |9j 
we present a slightly different version of this plot, where we 
have instead plotted the logarithm of the jet speed in order to 
alleviate crowding close to the x-axis. In calculating the lu- 
minosities, we have used the median total 2 cm VLBA flux 
density measured during the period 1994.0-2004.0 and have 
k-corrected the luminosities assuming a flat spectral index 
(a = 0). 

The dotted curve represents a jet with /i = 2.5 mas y" 1 , flux 
density of 1 .5 Jy, and flat spectral index, plotted for a range of 
redshifts. These values correspond to the maximum angular 
speed measurement and the flux density limit of the survey, 
respectively. The slope of this curve is roughly ~ 0.5 in this 
log-log plot, since the apparent velocity depends linearly on 
luminosity distance, whereas observed luminosity varies as 
luminosity distance squared (the k-correction causes a slight 
curvature at high z). As in the case of the maximum speed 
versus redshift plot (Fig. |8), the empty region between this 
curve and the data points indicates that the envelope is not the 
re sult of observation al measurement limits. 

ICohen et alj y007) pointed out that the shape of the enve- 
lope was very similar to the locus of points created by a single 
jet oriented over a range of different viewing angles. The free 
parameters involved in fitting this "aspect curve" to the en- 
velope are the jet's unbeamed luminosity (which determines 
its horizontal shift), its Lorentz factor (which determines its 
shape and vertical extent), and the Doppler boost index p, 
(which influences both the curve's shape and its slope at low 
apparent luminosities). 

For reasonably high values of the Lorentz factor, we can 
derive the approximate slope of the left hand side of the aspect 
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FIG. 8. — Plot of superluminal speed versus redshift for the fastest robust, non-accelerating component in 119 MOJAVE jets. Error bars smaller than the 
symbol sizes have been omitted. The red circles indicate quasars, the filled blue squares BL Lac objects, and the green triangles radio galaxies. The dashed curve 
represents an angular apparent speed of 2.5 mas y - ' . 



curve, where 2/T <; 5 « T (i.e., 9 » sin" 1 (1/T)): 

/3 a 2 pp = 2<5r-5 2 -l (3) 

~2(5r. (4) 

We also have 

L = dS p , (5) 

where C\ is a constant. Combining equations © and <j5j 
yields 

/3 app =V2T (L/Ci) (1/2 "\ (6) 

thus the slope in the log /3 app vs. log L plot is 0.5 p~ l . 

The solid curve in Figure [9] represents a jet with Lorentz 
factor T = 50 and an intrinsic unbeamed luminosity of 3 x 
10 25 WHz" 1 , viewed at angles ranging from 5 x 10" 4 to 20 
degrees, and Doppler boosted by a factor S 163 . The tick 
marks are located at one degree intervals moving from right 
to left, beginning at one degree from the line of sight. The 
aspect curve provides a reasonably good fit to the envelope, 
although there is considerable freedom in the fitting parame- 
ters. We determined the boosting index by performing a least 
squares fit to the ten data points between 10 24 7 W Hz" 1 < L < 
10 28 W Hz" 1 that lie along the upper envelope. This yielded a 
slope of 0.307 ± 0.015, which corresponds to p = 1 .63 ± 0.08. 



This Doppler boost is somewhat lower than that expected 
from a continuous flat-spectrum jet (p = 2 + a, where a = 0). 

Although the aspect curve appears to provide a good fit to 
the upper envelope, we find little justification for it in the 
framework of the relativistic beaming model. One interpre- 
tation would require that every jet located on the upper left 
edge of the envelope have a Lorentz factor of 50 and an iden- 
tical unbeamed luminosity. The lowest apparent luminosity 
jets in this sub-group would then have viewing angles of 10- 
20 degrees, and their Doppler factors would be on the order 
of unity. However, since we are dealing with a flux-density- 
limited beamed sample, there is a strong preference for high- 
Doppler factor jets. At a given observed luminosity, e.g., 
10 26 WHz" 1 , there is a much higher probability of having ajet 
with lower intrinsic luminosity and moderately high Doppler 
factor, in part because of the steepness of the parent lumi- 
nosity function, and also because of the rarity of ult ra-high 
Lorentz factor jets in the parent population (see Section l4.2.1b . 
Thus, there is no a-priori reason to expect a sharp upper en- 
velope bounded by a curve with T = r max . Another possible 
explan ation for the enve l ope, e xplored by iLister & Marscheri 
(119971) and ICohen et alJ (|2007), is that there is a correlation 
between intrinsic jet power and Lorentz factor, such that a jet 
with a particular T must have a certain minimum unbeamed 
luminosity L(T). 
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15 GHz VLBA Luminosity (W Hz ') 

FIG. 9. — Plot of maximum jet speed versus 2 cm VLBA luminosity for 1 19 jets in the MOJAVE survey. The red circles indicate quasars, the filled blue squares 
BL Lac objects, and the green triangles radio galaxies. Vertical error bars smaller than the symbol sizes have been omitted. The horizontal error bars indicate 
the range of luminosity of each source during the period 1994.0-2004.0, and the symbols correspond to the median values. The dotted curve corresponds to a jet 
with an apparent speed of 2.5 mas y - ' and a flux density of 1.5 Jy, plotted over a range of redshifts. The solid curve represents a jet at viewing angles ranging 
from 5 X 10 -4 degrees to 20 degrees, with Lorentz factor = 50 and an intrinsic unbeamed luminosity of 3 X 10 25 W Hz~' , Doppler boosted by a factor <5 163 . Ticks 
are plotted along the solid curve at one degree intervals, starting at a viewing angle of one degree and increasing towards the left. 
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FIG. 10. — Left panel: plot of apparent speed versus observed luminosity for 135 simulated jets (see Section 14.2.31 . The solid curve is identical to that in 
Figure [9] however, the plotted x-axis ranges of the figures are different. Right panel: identical simulation to the left-hand panel, but with a different random 
number seed. 
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We have attempted to investigate these possibilities by per- 
forming Monte Carlo simulations of a beamed parent popula- 
tion, from which we select randomly generated jets above the 
MOJAVE flux density threshold. The intrinsic jet luminosi- 
ties are generated assuming a power law pure-density evolu- 
tion LF with sl ope -2.65 and r edsh ifts between 0.04 and 4, 
as described by Cara & Lister (2008). Because of the small 
number of low-luminosity jets in the MOJAVE sample, these 
parameters are chosen to optimize the fit to the quasar pop- 
ulation. The jets are assigned random orientations, and their 
Lorentz factors are drawn from a power law distribution with 
slope = —1.5 and 3 < T < 50. The latter provides a reason- 
able fitto the overall distribution of maximum jet speeds (Sec- 
tion l4.2.lT ). Each jet is beamed using a Doppler boost index of 
p = 1.63. The simulation is provided a random number seed, 
and jets are generated until a flux-density-limited sample of 
comparable size to the MOJAVE survey is obtained. With 
the choice of parameters above, and 10 24 W Hz -1 < intrinsic < 
10 276 W Hz" 1 , this typically required ~ 20 million parent jets. 

In Figure [lOl we show the results of two simulations us- 
ing these parameters, albeit with different random number 
seeds. The fact that one of the simulations appears to repro- 
duce the envelope, while the other does not, reflects the statis- 
tical fluctuations that are inherent in selecting a biased subset 
of sources from a very large parent population. The appar- 
ent flux density of any particular jet is determined by the joint 
probability density of four variables, namely redshift, intrinsic 
luminosity, viewing angle, and Lorentz factor. Thus, there are 
multiple ways in which a jet can make it into a beamed sam- 
ple, such as having a low redshift, or a high Doppler factor, 
or both. Since the flux density limit of the MOJAVE survey is 
shallow, we are likely situated too far out in the statistical tail 
to fully constrain particular aspects of the parent population, 
such as an intrinsic relation between speed and luminosity. In 
order to address this issue, we are obtaining speed informa- 
tion on an extended sample of AGN jets that includes low- 
luminosity radio galaxies and new 7-ray blazars detected by 
the Fermi observatory. 

5. SUMMARY 

We have investigated the pc-scale kinematic properties of 
135 jets in the complete MOJAVE flux-density-limited survey 
of compact, radio-selected AGN. The observational dataset 
consists of 2424 VLBA images at 2 cm from our program and 
the VLBA archive between 1994 August 3 1 and 2007 Septem- 
ber 6, and represents a significant improvement over previous 
multi-epoch AGN surveys in terms of statistical complete- 
ness, temporal coverage, and overall data quality. 

Our main results can be summarized as follows: 

(i) By modelling individual bright jet features with Gaus- 
sian components, we were able to robustly identify and track 
a total of 526 moving features in 127 AGN jets. Three AGN 
(0235+164, 1324+224, 1741-038) had very little discernible 
mas-scale jet structure at any epoch, while five others dis- 
played too much scatter to reliably cross-identify their fea- 
tures across multiple epochs. We have also identified two 
bright AGN (0552+398 and 0642+449) with extremely sta- 
ble milliarcsecond-scale structures that have changed by less 
than 10 /ias y" 1 over a 12 year interval. 

(ii) We fit a simple linear model to the positions of each 
component over time in order to determine its speed and di- 
rection with respect to the (presumed stationary) core feature 
located near the base of the jet. Approximately 1/3 of the 



components had velocity vectors which did not extrapolate 
back to the core feature. A mo re detailed analysis presented 
elsewhere ( Ho man et alj 120091) links these 'non-radial' mo- 
tions to accelerations perpendicular to the velocity vector. 

(iii) The jet features show an overwhelming tendency to 
display outward motions, away from the core (98.5% of all 
components). Our analysis revealed only 8 inward-moving 
components, and in only one case (ID=2 in 2351+456) is the 
motion highly superluminal. VLBA data obtained since 2007 
on the latter feature has indicated that its inward motion has 
stopped, suggesting that it may be a case where the jet has 
crossed our line of sight. 

(iv) Features within an individual jet tend to have a mod- 
estly narrow distribution of apparent speed (three times 
smaller than the overall speed dispersion in the sample), sup- 
porting the idea that they are connected with the intrinsic 
flow. We identified only 3 1 features whose motions were slow 
and significantly smaller than others in the same jet. These 
"low-pattern speed" features are more prevalent in the BL 
Lac objects (10 of 22) than in the quasars (15 of 86), and 
are rarely found at de-projected distances greater than 50 pc 
downstream of the core. They may represent strong, stable 
re-collimation shocks of the type commonly seen in numeri- 
cal jet simulations. 

(v) We have analyzed the kinematics of five MOJAVE 
sources that have jets of similar brightness located on either 
side of the central AGN. We find that at least two (NGC 1052 
and Cygnus A) of these two-sided jets display a close simi- 
larity in the ejection epochs and apparent speeds of their jet 
and counter] et, providing additional evidence that the moving 
features are related to the flow. 

(vi) Nearly 60% of the 526 robust components in the sample 
had at least 10 epochs of measurement; for these we also per- 
formed a two-dimensional fit to determine their velocities and 
accelerations. Approximately one third of these display sig- 
nificant accelerations at the 3er level, indicating that changes 
in speed and/or direction are a common feature of AGN jets. 

(vii) We have compiled a maximum speed statistic for each 
jet, which is equal to the fastest robust feature observed in the 
source. The distribution of these maximum speeds is peaked 
at /3 a pp — 10, and falls off sharply toward higher values. This 
is indicative of a parent population that has a distribution of 
bulk Lorentz factors that is not peaked at high values. The 
fastest jet in our sample (0805-077) contains a component 
moving at 50.6 ±2.1 c, which is likely representative of the 
upper end of the AGN jet Lorentz factor distribution. 

(viii) The maximum jet speed increases with redshift out to 
Z—l, which reflects a selection bias towards higher Lorentz 
factor jets at higher redshifts. At closer distances, the jets do 
not require as high beaming factors to meet the flux density se- 
lection criteria, and the probability finding a high luminosity, 
high Lorentz factor jet is much lower because of the smaller 
co-moving volume available. 

(ix) We have re-examined the upper envelope i n the /3 app vs. 
appare nt luminosity plane first reported by ICohen et alj 
(2007) using our new data set. The envelope is not the re- 
sult of any observational measurement bias, and is very sim- 
ilar in shape to an "aspect curve" on the plane traced by a 
single high-Lorentz factor jet as it is rotated through differ- 
ent viewing angles. The statistical probability of having this 
envelope, however, is low, in terms of the standard beaming 
model, which we have confirmed using Monte Carlo simula- 
tions with multiple random number seeds. It is possible in- 
stead that the envelope arises because of an intrinsic relation 
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between jet speed and luminosity in the parent population, 
although deeper samples will likely be needed to fully inves- 
tigate this scenario. 
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Jet Kinematics 
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NOTE. — Columns are as follows: (1) IAU name (B1950.0); (2) component number; (3) number of fitted epochs; (4) mean flux density; (5) mean distance from core component; (6) mean position angle 
with respect to the core component; (7) magnitude of fitted angular velocity vector; (8) fitted speed in units of the speed of light; (9) position angle of velocity vector; (10) offset between mean position angle 
and velocity vector position angle; (11) acceleration parallel to velocity direction; (12) acceleration perpendicular to velocity direction; (13) fitted ejection date; (14) date of reference (middle) epoch used for 
fit; (15) right ascension error of individual epoch positions (mas); (16) declination error of individual epoch positions (mas). 
a Component shows significant acceleration. 

Component shows significant non-radial motion. 
c Component shows significant inward motion. 
^ Component is unusually slow (low pattern speed) (see Section [3~5] , 



